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' We discuss several open problems of fusion reactions induced by weakly bound nuclei. For 

I this purpose, we solve a one dimensional three-body Hamiltonian with the coupled-channels 

, formalism. We show that the continuum-continuum couplings substantially reduce the total 

K ■ fusion probability at energies above the barrier compared with the no-breakup case, although 

the fusion probability remains enhanced at subbarrier energies. We then discuss a role of 
transfer process in fusion of weakly bound nuclei, and point out that removing spurious 
0^ . Pauli forbidden transfer components from the calculation may be crucial at energies below 

the barrier. Calculations based on the three-body classical trajectory Monte Carlo (CTMC) 
method are also presented in order to discuss how to model complete fusion process. 
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§1. Introduction 



I The recent availability of radioactive beams has opened up the opportunity to 

' study the interactions and structure of exotic nuclei close to the drip lines. Unstable 

neutron-rich nuclei having very weakly bound neutrons exhibit characteristic features 
such as a neutron halo and a low energy threshold for breakup. These features may 
dramatically affect fusion and other reaction processes. 
O , One of the main issues in the study of reaction induced by weakly-bound nuclei is 

^ I to clarify how the breakup process influences fusion reactions at energies close to the 

Coulomb barrier. ^^^^''^ In Ref. 4), we have performed coupled-channels calculations 
in order to address this question by discretizing in energy the particle continuum 
^ ' states of the projectile nucleus. Defining the complete fusion as absorption from 

■ bound state channels, we have demonstrated that fusion cross sections are deter- 

mined by the competition of two mechanisms, i.e., dynamical modulation of fusion 
barrier and flux loss due to the breakup couplings. Their net effect differs depending 
on the bombarding energy: the complete fusion cross sections are enhanced at en- 
ergies below the barrier, while they are hindered above the barrier, compared with 
cross sections for a tightly bound system. On the other hand, the total fusion cross 
sections (a sum of complete and incomplete fusion cross sections) are enhanced at 
subbarrier energies while they are close to the no-breakup calculations at above. In 
these calculations, we did not include the couplings among the continuum channels. 

In this contribution, we would like to discuss how these conclusions may be al- 
tered by more refined calculations. We particularly pay attention to the definition of 
complete fusion, the role of continuum-continuum couplings, and the effect of trans- 
fer channels. To this end, we solve a simple one dimensional three-body Hamiltonian 
with the coupled-channels method. As we will show below, our main conclusions for 
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the complete fusion remain valid at least in a qualitative way, although we find that 
the total fusion cross sections are hindered at energies above the barrier due to the 
continuum-continuum couplings. 

§2. Separation of complete and incomplete fusion 

Theoretically, the complete fusion 
refers to the capture of the whole projec- 
tile, while the incomplete fusion is a pro- 
cess where only a part of the projectile is 
captured. In our previous calculations,^^ 
we defined complete fusion as absorp- 
tion from bound state channels, while 
the incomplete fusion as those from con- 
tinuum state channels in the projec- 
tile. The basic idea was that a contin- 
uum state is a doorway of breakup, and 
thus is related to incomplete fusion pro- 
cess. However, when the breakup takes 
place inside the Coulomb barrier, all 
the breakup fragments would likely be 
captured by the target nucleus, leading 
to complete fusion reaction. This 
breakup followed by complete fusion was 
not taken into account in the definition 
of complete fusion in Rcf. 4), and the 
previous calculation provides only the 
lower limit of complete fusion cross sec- 
tions. 

In order to estimate the contribution of the breakup followed by complete fusion, 
one of us has developed a three body classical trajectory Monte Carlo (CTMC) 
method.^-* For a three body Hamiltonian that consists of the target (T) and two 
projectile fragments (Pi and P2), two-dimensional classical Newtonian equations are 
solved to obtain the time evohition of the coordinates and velocities of the fragments. 
The initial conditions are that the projectile, with its two fragments in random 
orientation, starts far from the target with impact parameter b. As the projectile 
moves towards the target, three processes are possible, depending the value of the 
impact parameter b: (i) the projectile as a whole or both of the fragments are 
absorbed by the target, (ii) only one fragment is absorbed, and (iii) neither fragment 
is captured. We assume that the breakup occurs when the distance between the 
projectile fragments exceeds their potential barrier radius, and that the projectile 
fragment Pi is absorbed by the target nucleus when the relative distance between 
the target and the fragment is smaller than Tabs = 1-1 x {A^^ + ^p^^)- See also Ref. 
11) for details. 
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Fig. 1. Comparison of the thrcc-body classical 
trajectory Monte Carlo calculations with 
the experimental data for ®Li -t-^^^Bi reac- 
tion. The solid and the dotted lines arc 
complete and total fusion cross sections, 
respectively. The contribution of the ab- 
sorption of the whole projectile from bound 
states is denoted separately by the dashed 
line. The experimental data are taken from 
Ref.«) 
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Figure 1 shows the results of the three body CTMC method for the fusion 
reaction of ^Li +^°^Bi.^) The dotted hne denotes the sum of complete and incomplete 
fusion cross sections. Since this is a classical calculation, cross sections are finite only 
at energies above the Coulomb barrier. The dashed line indicates the absorption of 
the whole projectile from the bound states, while the solid line shows the total 
complete fusion cross sections. The difference between the solid and the dashed line 
thus represents the contribution from the breakup followed by the complete fusion 
process. As we see, there is a significant contribution of this process at energies 
above the barrier. This process has not been considered in any coupled-channels 
calculations so far, and the present result suggests that a consistent definition of 
complete fusion is necessary when one compares experimental data with theoretical 
calculations. 

§3. Role of continuum-continuum couplings 

Let us now discuss the effect of continuum-continuum couplings in the coupled- 
channels calculations. Diaz-Torres and Thompson showed that this effect increases 
the irreversibility of breakup process, and thus hinders the total fusion cross sec- 
tions. In fact, their calculation suggests that the total fusion cross sections are 
smaller than the no-coupling case at energies above the barrier, although those are 
larger at energies below the barrier. A similar hindrance has been found in the time 
dependent wave packet approach of Yabana et al.^^'"^^ Their calculations indicate 
that the total fusion probability is smaller than that for the core-target system, sug- 
gesting that it is much smaller than the prediction of the folding model calculation, 
which statically takes into account the halo structure of the projectile nucleus. They 
attribute this hindrance of fusion probability to a spectator role of the valence neu- 
tron. In this section, we revisit this problem and discuss whether the hindrance 
of total fusion probability observed in the time dependent wave packet approach 
can be understood in terms of the effect of continuum-continuum couplings in the 
coupled-channels approach. 

To this end, we use a one dimensional three-body Hamiltonian which was used 
by Yabana and Suzuki for their time-dependent wave packet calculation for fusion of 
^°Li -I- ^°Ca system.^) Assuming that the projectile consists of a core nucleus (C) and 
a valence neutron (n), and that the target mass is infinite, the model Hamiltonian 
reads^^ 

where X and x are the relative coordinate between the target and the center of mass 
of the projectile, and that between the core and the valence neutron in the projectile, 
respectively. AI = rrin + mc is the total mass of the projectile, while /x = nin-mc/M 
is the reduced mass for core-neutron system. We use the same potentials with the 
same values for the parameters as in Ref.6), except for the Xm parameter in the 
Coulomb interaction which we somewhat increase in order to increase the barrier 
height. 
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In order to solve this Hamilto- 
nian with the coupled-channels frame- 
work, we first obtain the eigenstates 
for the projectile Hamiltonian, i^pioj = 
+ ^ncix), by expanding the 
wave functions on the orthogonal poly- 
nomial basis, which was recently advo- 
cated by Perez-Bernal et al}"^^ This ba- 
sis is constructed from the ground state 
wave function of the Hamiltonian 0o(x) 
by multiplying a polynomial function to 
it {<t)Q,x(l)Q,x^(f)Q^ - ■ ■ ,x"'(^o)- Since the 
wave function is expanded on the dis- 
crete basis, the continuum states, ob- 
tained by diagonalizing the Hamilto- 
nian, are also discrctizcd in energy. We 
use the polynomial basis up to n=28. 
The energy for the first five states are 
given by: eo = -26.09 MeV, ei = 
-11.42 McV, €2 = -0.515 McV, eg = 
0.135 MeV, and £4 = 0.193 MeV. Fol- 
lowing Ref. 6), we assume that the 
ground and the first excited states in 
this Hamiltonian are occupied by the 
neutrons in the core nucleus, and con- 
sider the second excited state as the 
ground state of the core-neutron system, 
4>gs{x). 

Figure 2 shows the total barrier 
penetrability obtained with the coupled- 
channels method (the solid line) in the 
linear (the upper panel) and in the log- 
arithmic scale (the lower panel). We in- 
clude the first seven discretized continuum states together with the ground state. At 
low energies, some of the continuum states arc kincmatically forbidden (i.e., closed 
channels), and it is not easy to obtain numerically stable solutions of the coupled- 
channels equations. We avoid this problem by using the modified log-derivative 
method. This method is in fact very powerful, and we have confirmed that our re- 
sults are stable against the mesh size and the matching radius. We have also checked 
that the low energy penetrability does not change much even when we include more 
continuum state channels in the coupled-channels equations. The penetrability still 
shows a resonance-like structure. This occurs at a threshold energy for each dis- 
cretized channel, and can thus be attributed to the threshold effect. 

The result of the coupled-channels calculations is compared with the fusion prob- 
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Fig. 2. The total penetrability for the one 
dimensional three body Hamiltonian ob- 
tained by the several methods. The solid 
and the dashed lines are the results of 
the coupled-channels calculations with and 
without the continuum-continuum cou- 
plings, respectively. The dot-dashed line is 
the penetrability for the folding potential, 
while the dotted line denotes the penetrar- 
bility for the core-taget system. 
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ability (barrier penetrability) for the folding potential (the dot-dashed line), 
^foid(^) = J dx(f>gs{xf 

We find that the total penetrability is significantly hindered by the continuum cou- 
plings at energies above the barrier height of the folding potential, although it is still 
enhanced compared with the penetrability for the folding potential at energies below 
the barrier (see the lower panel). Notice that the total penetrability is even smaller 
than the penetrability for the core-target system (the dotted line) at energies above 
the barrier. This is reminiscent of the result of the time dependent wave packet 
approach.^) 

The result without the continuum-continuum couplings is denoted by the dashed 
line. As was pointed out in Ref. 5), this calculation significantly overestimates 
the barrier penetrability compared with the full coupled-channels calculation. A 
part of the reason why the continuum-continuum couplings reduce the penetrability 
can be understood as follows: the wave function of the (discrctizcd) continuum 
state considerably extends outside the barrier. The diagonal matrix element of the 
coupling Hamiltonian for a continuum state channel is thus much smaller than that 
for the ground state, eq. (3-2). Notice that the diagonal matrix element is given by 
the folding potential itself for all the channels if the continuum-continuum couplings 
are neglected. The change of the diagonal matrix elements can be regarded as 
"effective Q- value", which is large and positive around the barrier region. If the 
coupling is strong enough, the penetrability can then be reduced due to the so called 
anti-adiabatic effect. ^^^'^^^ 

§4. Role of transfer process 

We next discuss the effect of transfer couplings on fusion of a halo nucleus. 
Compared with the breakup effects, a discussion on the transfer effect has yet been 
scarce. The time dependent wave packet approach of Yabana et al. shows that 
the fusion probability is hindered over that for the core-target system even at en- 
ergies well below the barrier. This conclusion is qualitatively different from the 
result of the coupled-channels approach,^)' ^) which predicts an enhancement of fu- 
sion probability at low energies, as was discussed in the previous section. In both 
of these approaches, the multipole decomposition is usually introduced to the total 
wave function. Yabana et al. examined the effect of the angular momentum trunca- 
tion for the neutron-core system, and showed that the fusion probability is enhanced 
when the angular momentum is limited to low values.^) A high value for the angu- 
lar momentum is required in order to simulate the transfer process of the neutron 
to the target nucleus, and Yabana et al. claimed that the enhancement of fusion 
cross sections found in the coupled-channels approach is an artifact of neglecting 
the transfer coupling. They needed up to Z = 70 in order to get a converged result. 
Diaz- Torres and Thompson also studied the effect of angular momentum truncation 
in the coupled-channels framework.^) Their result seems to indicate that the fusion 
cross sections indeed get smaller when higher anuglar momentum components are 
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included in the calculation (sec Fig. 6 in Ref. 5)). 

However, we here ask a general 
question: is it really physical to include 
all the partial wave components for the 
neutron-core system? It is true that 
one can in principle solve the three-body 
Hamiltonian exactly if one includes all 
the partial waves. In such a calculation, 
all the transfer channels are automati- 
cally entered. Notice that these include 
even a transfer to a deep bound state of 
the target nucleus, which must be for- 
bidden by the Pauli principle. Those 
Pauli forbidden transfer channels most 
likely have a large positive Q value, and 
the fusion cross sections may be artifi- 
cially hindered due to the anti-adiabatic 
effect^^^ if these channels are included 
in the calculation. Also, when the ini- 
tial configuration of the valence neutron 
is at an excited state of the Hamiltonian 
due to the Pauli principle, the forbidden 
dc-cxcitations to lower states in the pro- 
jectile may be entered unless those tran- 
sitions are explicitly excluded by e.g. 
the basis expansion. Those forbidden 
de-excitations may also lead to an ar- 
tificial reduction of fusion cross sections 
for the same reason. 

In order to examine the effect of the 
spurious Pauli forbidden transfer chan- 
nels in the projectile, we again solve the 
one dimensional three body model in 
the previous section with the coupled- 
channels method by including several 
transfer channels. For the parameter set 
specified in the previous section, there 
are five bound states in the target nu- 
cleus at -40.5, -31.9, -20.6, -8.85 and 
-0.357 MeV. We assume that the 
ground state as well as the first and the 

second excited states in this Hamiltonian are already occupied. The solid line in 
the upper panel of Fig. 3 is the result of the coupled-channels calculation where 
only the transfer channels to the two allowed states in the target are included. We 
compute the transfer coupling matrix elements in the no-recoil approximation using 
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Fig. 3. The effect of the Pauh forbidden trans- 
fer channels on the transmission (the up- 
per panel) and the reflection (the lower 
panel) probabilities. The dashed and the 
solid lines in the upper panel are obtained 
by solving the coupled-channels equations 
with and without the Pauli forbidden chan- 
nels, respectively. The dotted line in 
the upper panel is the penetrability for 
the folding potential. The solid lino in 
the lower panel is the reflection probabil- 
ity for the incident channel obtained by 
the coupled-channels calculation with the 
Pauli forbidden channels. The dashed and 
the dot-dashed lines are the transfer proba- 
bilities for the highest and the second high- 
est bound states in the target nucleus. The 
total reflection probability for the Pauli 
forbidden channels is denoted by the dot- 
ted line. 
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the neutron-target potential VnT as the transfer interaction. In order to mimic the 
recoil effect, we multiply a factor of 2{mc + l)/{2mc + 1) to the transfer matrix 
elements. ^'^^'^^^ For simplicity, the continuum excited states in the projectile are 
not included. In this calculation, one can clearly see that the fusion probability is 
enhanced due to the transfer couplings at energies below the barrier compared with 
the folding model calculation (the dotted line) . The result of the coupled-channels 
calculation with the Pauli forbidden channels is given by the dashed line. One sees 
that the enhancement of fusion probability is much less at energies below the barrier 
due to the anti-adiabatic effect. 

The lower panel of Fig. 3 shows the reflection probabilities for the elastic (the 
solid line) and the transfer channels. The dashed and the dot-dashed lines repre- 
sent the transfer probability for the highest and the second highest bound states 
in the target nucleus, respectively. Due to the Q-value matching, the total transfer 
probability is dominated by the former, as has been found in the time dependent 
wave packet approach.^) The dotted line is the total reflection probability for the 
Pauli forbidden channels. Notice that this probability is much smaller than the 
transfer probability for the allowed channels, and may be negligible compared to the 
latter. Yet, the Pauli forbidden processes influence the barrier penetrability (the 
fusion probability) in a significant way. We thus conclude that care must be taken 
in discussing the angular momentum truncation, especially at energies around the 
barrier. 

§5. Summary and discussion 

Two alternative approaches have been used in the literature in order to discuss 
the subbarrier fusion reaction induced by a weakly bound nuclei. A traditional 
approach is to solve the coupled-channels equations by discretizing in energy the 
continuum states in the projectile nucleus (CDCC). The second approach employs 
the time dependent dynamics and solves the time evolution of a wave packet. If 
one neglects the dynamical effects, the total fusion probabilities would have been 
larger than that for the core-target system, since the barrier height is considerably 
reduced due to a large extension of neutron wave function in the projectile nucleus. 
On the contrary, the time dependent approach predicts smaller fusion probabilities, 
suggesting that there is a huge dynamical effect which compensates the static effect. 
Using a simple one dimensional three body Hamiltonian, we have demonstrated that 
this result can be understood in terms of the effect of couplings among the continuum 
states. In this conference, Signorini et al. reported that the measured fusion cross 
sections for the -'^^Be+^'^^Bi are similar to those for the ^'^Be+^'^^Bi system despite 
the halo structure of the ^^Be nucleus. It might be that the enhancement of fusion 
cross sections due to the static effect (reduction of the barrier height due to the halo 
structure) is compensated by the dynamical breakup effect, resulting in similar cross 
sections between the two systems. 

At energies below the barrier, the two theoretical approaches predict a contra- 
dictory result to each other. In the coupled-channels approach, the fusion probability 
is found to be larger than the prediction of the folding model. In contrast, the time 
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dependent wave packet approach predicts a hindrance of fusion probabihty even at 
energies below the barrier. We have argued that care must be taken in the partial 
wave decomposition for the wave functions because spurious Pauli forbidden transfer 
channels are inevitably coupled when all the partial waves are included. We have 
found that the Pauli forbidden transfer couplings significantly hinder the total fusion 
probability even if the transfer probability itself may be negligible. 

In order to draw a more definite conclusion on fusion cross sections for a weaJily 
bound system at energies below the barrier, further calculations would be required 
by including only the Pauli allowed transfer channels together with the breakup 
channels. The time dependent wave packet approach would need to introduce the 
projection operator in order to exclude explicitly the Pauli forbidden channels. In 
the coupled-channels approach, the angular momentum between the core and the 
valence neutron should be limited to low values in order to avoid the Pauli for- 
bidden transfer, and should explicitly include the transfer couplings in addition to 
the breakup channels. The combined effect of breakup and transfer in the coupled- 
channels approach has not been studied so far in the context of subbarrier fusion, 
and it would be an interesting future work. 

A further complication arises from a separation of complete and incomplete 
fusion cross sections. Using the three body classical trajectory Monte Carlo method, 
we have demonstrated that there is a significant contribution of the breakup followed 
by complete fusion process to the total complete fusion cross sections. It would be 
another interesting future problem to model this process using a quantum mechanical 
model for subbarrier fusion of unstable nuclei. 
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